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The generic bound of quantum speed limit time (the minimal evolution time) for a qubit system interacting 
with structural environment is investigated. We define a new bound for the quantum speed limit. It is shown that 
the non-Markovianity and the population of the excited state can fail to signal the quantum evolution accelera¬ 
tion, but the initial-state dependence is an important factor. In particular, we find that different quantum speed 
limits could produce contradictory predictions on the quantum evolution acceleration. 
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I. INTRODUCTION 

The maximal dynamical speed of quantum system is a fun¬ 
damental concept in many areas of quantum physics, such as 
quantum communication [1,2], quantum metrology [3,4], op¬ 
timal control [5], etc. The quantum speed limit (QSL), defined 
as the minimal evolution time between two states was first in¬ 
troduced by Anandan-Aharonov using the Fubini-Study met¬ 
ric [6]. Later, the unified QSL is given for closed systems 
subject to the unitary evolution, by the combining the con¬ 
clusion of the variance of the energy, i.e., Mandelstam-Tamm 
bound (MT-QSL) [6-10] and the average of the energy, i.e., 
Margolus-Levitin bound (ML-QSL) [11-13], which is rig¬ 
orously described by = max{7r/i/(2A£)), 7r?i/(2£’)}. 

The QSL bound for the driven system [14-17], and for the 
mixed quantum state evolution [18-20] were also investigated. 
In reality, the system interacts with the its environment in¬ 
evitably, so the theories of open quantum systems are usually 
employed [21]. Recently, the QSL for open system (nonuni¬ 
tary evolution) has attracted intensive interests, such as char¬ 
acterizing QSL for open system using quantum Fisher infor¬ 
mation [22], relative purity [23], Bures angle [24], and so on 
[20, 25-29]. The works [22, 23] also apply to non-Markovian 
dynamics. In particular, in Ref. [24], the authors extended the 
ML-QSL from closed system to open system using the tech¬ 
nology of the von Neuamnn trace inequality [30, 31] and op¬ 
erator norm [32], and showed that the non-Markovian effects 
[33-37] can speed up the quantum evolution and the ML-QSL 
is tight for the open system. In Ref. [28], the authors showed 
that the mechanism of the acceleration in open systems is de¬ 
termined not only by the non-Markovianity but also by the 
population of excited states during the quantum evolution. 

In this paper, we find that the initial-state dependence is also 
a key factor to the acceleration of quantum evolution. In par¬ 
ticular, the predications for the acceleration of quantum evolu¬ 
tion based on different QSLs could be quite inconsistent. Here 
we first define a tight bound for the quantum speed limit with¬ 
out using the von Neumann trace inequality (We denote it by 
NI-QSL). Then we use the ML-QSL and the NI-QSL to inves¬ 
tigate the damped Jaynes-Cumming (J-C) model with a super¬ 


position state as the initial state and the dephasing model. It 
is found that neither the non-Markovianity nor the population 
of the initial excited state is competent for the acceleration of 
quantum evolution, but the different initial states bring great 
effects. In addition, the ML-QSL and the NI-QSL sometimes 
produce contradictory predictions on the evolution accelera¬ 
tion, even though they are consistent in some particular cases. 
This implies that the QSL could need us further consideration 
and investigation. This paper is organized as follows. In the 
Sec. II, we give a brief introduction about the ML-QSL and 
define the NI-QSL. In Sec. Ill and IV, we study the damped 
J-C model and the dephasing model, respectively. Some dis¬ 
cussions and conclusion are drawn at the end. 


II. THE QUANTUM SPEED LIMIT 

In Ref. [24], the authors derived the ML-type quantum 
speed limit bound for a pure state |'0o) based on the von Neu¬ 
mann trace inequality and operator norm in the open system 
through geometric approach. The ‘distance’ between the ini¬ 
tial state po = IV'o)('0o| and its target state pt is measured by 
Bures angle, which reads 

^Po,Pt) = arccos(y/(^/>o|pt IV'o))- (1) 

The generalized time-dependent nonunitary equation can be 
expressed as pt = Lt (pt)- The time derivative on the Bures 
angle leads to the relation 

2cos (£) sin (£) £ < |('0o| Pt IV'o)! = (pt) Poll • (2) 

Based on the von Neumann trace inequality, the right-hand 
side of Eq. (2) can arrive at 

n 

\tr{Lt (pt) Poll < X (3) 

i=l 

with (Ti i and <72,i being the singular values of matrix Lt (pt) 
and po in decreasing order. Integrating Eq. (2) over time, one 
can obtain the ML-QSL for the nonunitary generator as 

ml ^ sin^[£(po,pt)] 

^QSL ^ ^op 
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( 4 ) 
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where = ( 1 /t) J* dr || At(Pr)|lop with IHI being the op¬ 
erator norm (the maximum singular value) of the matrix. 

However, if the von Neumann trace inequality Eq. (3) and 
operator norm are not used, we can directly obtain, through 
integrating Eq. (2) over time, 

-, (5) 

(lA)/o dT\{^ q\L t{pr) |^/>o)| 

which we call the NI-QSL to distinguish from the ML-QSL. 
The prominent feature of Eq. (5) is that it provides a tighter 
bound than ML-QSL given by Eq. (4). Our main conclusion 
will be obtained by comparing the two types of QSL for the 
dynamics of open quantum systems. 


III. THE QUANTUM SPEED LIMIT TIME EOR THE 
DAMPED JAYNES-CUMMING MODEL 


The first model we will consider includes a single qubit in¬ 
teracting with vacuum reservoir [21], which is also called the 
damped J-C model. The whole Hamiltonian of the system 
and reservoir is iJ = uJoa^ + J^k^kblbk + J2ki9kcr+bk + 
gl.a-b\). The dynamics of the system can be described by 

Lt{pt) = Jt{cr-P(T+ - ^cr+pcr_ - ^a-pa+). (6) 

The structure of environment is assumed as the Lorentzian 
form: 

= Imp S(u„ - u,) = ^ ^ y ■ (7) 

where A is the spectral width of the reservoir and 70 is the 
decay of the system. The Markovian and non-Markovian 
regimes can be distinguished by the relation of parameters 7 
and A. In the Markovian regime, we have 70 < A/2 and in the 
non-Markovian regime, we have 70 > A/2 [33-37]. 

If the initial state of system is a superposition state 

I'!/)) = ae*^|I) -I- v/l — a^jO), (8) 

the final reduced state after the evolution will be [21] 

„ _ f aVl - \ 


where q{t) is determined by the integro-differential equation 
q{t) = Jp dTf{t — T)q{T) with the correlation kernel re¬ 
lated to the spectral density of the reservoir as /(f — r) = 
J Using the Laplace transformation 

and its inverse transformation, q(t) can be given by 


q{t) = e 2 



^ • u 

+ — Sinn 

a 



( 10 ) 


with d = The time dependent decay rate 74 

in the Eq. ( 6 ) is given by 74 = —Im ( ) . After some 


calculations, we can obtain the ML-QSL for the initial state 

|' 0 ) as 


,ML _ 
^QSL - 


f|a|(l-g(Q)[l-(l-2a2)g(f)] 


fo \/l - (1 - 4g(T)2)a2g(r) 


dr 


( 11 ) 


Erom Eq. (11), one can find that the ML-QSL is related 
to the a, q(t) and q{t). The variation of ML-QSL with 
the parameter 70 and a is plotted in Panel (a) of Eig. 1. In Eig. 
1, the actual evolution time is chosen as f = 1 and A = 15. 
tqgL is symmetry about the population parameter a = 0 , so 
we just study the part a > 0. It is obviously shown that for 
7 o < A/2, the quantum speed limit time is below the actual 
evolution time 1. Namely, the quantum evolution displays the 
acceleration in the Markovian regime. This is quite different 
from that the acceleration only appears in the non-Markovian 
regime due to non-Markovianity which was discussed a = 1 
in Ref. [24]. One could imagine that the non-Markovianity is 
not the only reason. The population of the excited state is also 
a factor for the such an acceleration as that given in Ref. [28], 
which is determined by 

t 

where N {■) is the degree of non-Markovianity for dynamics 
defined as the total backflow of information [34] and |( 7 (f)| 



FIG. 1: (Color online). The quantum speed limit time for the initial 
state given by Eq. (8). Panel (a) is the ML-QSL and Panel 

(b) is the NI-QSL igsi,, where the parameter A = 15 and the actual 
evolution time is f = 1. 
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FIG. 2: (Color online). The quantum speed limit time for the initial 
state \^p) = (|1) + |0))/'\/2. Here A = 15. The dashed line is 
the actual evolution time f = 1 , the dash-dotted line stands for the 
ML-QSL and the solid line represents the NI-QSL fgsi,. 


is the population of the excited state, one should notice that 
in the Markovian regime, N (•) = 0, it will directly elimi¬ 
nate the role of the population \q{t)\^ in Eq. (12). Thus one 
possible reason for the acceleration subject to the ML-QSL is 
the dependence of the initial state. The details on the non- 
Markovianity is given by the Appendix A. 

If we use the NI-QSL to characterize the quantum speed 
limit, we have to calculate the NI-QSL for the initial state \ijj) 
as 


fo |2[1 - - (1 - 2g(r)2)a2]q(r)| dr 

The variation of ML-QSL with the parameter 70 and 
a is plotted in Panel (b) of Fig. 1. One can find that the 
quantum speed limit time is fixed to the actual evolution time 
in the Markovian regime, while the evolution acceleration is 
demonstrated in the non-Markovian regime, i.e., 70 > A/2. 
In order to give a more intuitive illustration, the QSLs 
and fgsi for the system with initial state |i/;) = ^(| 1 ) -|- | 0 )) 
are plotted in Fig. 2. In Fig. 2, the parameters are chosen 
the same as Fig. 1. The initial-state dependence can also been 
found from Fig. 1 by different a. 

Comparing the quantum speed limit and one 

can find that for the initial state 11 ) (i.e., a = 1 ). 

However, if a 1, and tqgi^ usually demonstrates dif¬ 
ferent behaviors which has been analyzed previously. Thus 
the ML-QSL and the NI-QSL produce contradictory predic¬ 
tions on the evolution acceleration. 


IV. QUANTUM SPEED LIMIT EOR THE DEPHASING 
MODEL 


In the following, we will consider another exactly solvable 
model, a two-level system coupling with a harmonic oscillator 
reservoir, it is also called the dephasing model [21]. In the 
Schrodinger picture, the total Hamiltonian is taken to be iT = 
+Y.k^kh\bk + Y.k^^{9kb\ + glbk)- The evolution 


operator of the system is Lt{pt) = ^{o'zPtO'z — Pt)- The 
initial state is 

!</) =^e*®|l) + ^r^|0), (14) 


and the dynamics of the reduced system is expressed as [ 21 ] 


Pt 


(3^ \ 


Taking the continuum limit of the bath mode and introduc¬ 
ing the spectrum J{uj) of the environment, we can find the 
dephasing factor 7 (f) given by 

7 Q = / dujJ{uj)coth{-—) - 2 -■ (16) 

Jo 2kBT 

The spectrum of the environment is chosen as the Ohmic-like 
spectrum with soft cutoff [38, 39] 

T(w) = expj-w/wc), (17) 

LJc 

where Wc is the cutoff frequency, rj is the dimensionless cou¬ 
pling constant and the parameter s > 0. For simplicity, we 
will assume that the cutoff frequency Wc is 1 . s determines the 
property of the environment such as the sub-Ohmic reservoir 
for s < 1, the Ohmic reservoir for s = 1 and the super-Ohmic 
reservoir for s > 1. Under the condition T = 0, f > 0 and 
s > 0, the dephasing factor 7 (f) can be obtained [40] by 


7(f) = 77 


cos[(s — 1 ) arctan(f)]r(s — 1 ) 

(l-pf2)(s-l)/2 


(18) 


where r(-) is the Euler Gamma Function. Thus the ML-QSL 
for the state in Eq. (14) can be given by 


.ML _ 
^QSL — 


2/3^1-/32(l-e-T'l*))f 
Jo |e“T'('^)7'(r)I dr 


(19) 


From Eq. (19), one can find that the ML-QSL tqgj^ de¬ 
pends on the dephasing rate 7 '(t), which is the time derivative 
of the dephasing factor 7 (f) and determined by 

y(i) = / (20) 

J w 

Here, we have used the zero temperature condition. In addi¬ 
tion, one can see that, if ^'{t) > 0 , the integral in denominator 
of Eq. (19) can be given analytically as 1 — It will ob¬ 

viously lead to that tqgj^ = 2/?-\/l — f3H. In other words, the 
quantum speed limit time is fixed for the given initial state. It 
is irrelevant to the property of the environment. However, in 
the region of 7 '(f) < 0, in usual the denominator of Eq. (19) 
is hard to give explicitly. So it is important to tell the sign of 
7 '(f) in order to give the ML-QSL. For the Ohmic-type spec¬ 
tral, the dephasing rate 7 '(t) can be given analytically as 

j'(t) = ry(l -I- sin(s arctan(f)). (21) 



















One can find that 7 '(f) is always positive if sin(s arctan(f)) > 
0 , i.e., s and t satisfy 
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k-K < sarctan(f) < {2k + l) 7 r, (22) 

otherwise, j'{t) is negative. A vivid illustration of Eq. (21) is 
provided in Fig. 3. The green region corresponds to the pos¬ 
itive 7 '(t), and the other black region stands for the negative 
7'(t). 

The ML-QSL versus s > 2 and p is plotted in Fig. 

4, where the lower layer corresponds to 2/3i/l — = 0.5 

and the upper layer corresponds to /3 = It can be found 

that the MF-QSF for /3 = -^ serves the tightest MF-QSF for 
all possible (3. Here we choose the actual evolution time as 
t = 3. Since arctan(f) € [0, 7 r/ 2 ), if s < 2, 7 '(t) will always 
be positive, which can be seen from Fig. 3. In this case, the 
MF-QSF will be a constant and independent of s and p, for 
example, = 1.5 for the lower layer and = 3 is 
just the actual evolution time for the upper layer. So the part 
for s < 2 is not shown in Fig. 4. Compared with Fig. 3, one 
can find that the MF-SQF keeps invariant until s ~ 2.5. For 
all s the MF-SQF given by the lower layer shows the acceler¬ 
ation of the quantum evolution, but no acceleration is shown 
by the upper layer for s < 2.5 . In both layers, the popula¬ 
tion of the excited state is not changed due to the dephasing 
dynamics. Therefore, the acceleration with different degrees 
could be understood by the non-Markovianity. In Ref. [39], 
it is shown that for s > 2 the dynamics for t ^ oo enter the 
non-Markovian regime. In fact, if t is finite, one can easily 
find that the negative 7 '(f) will lead to the non-Markovian dy¬ 
namics (see Appendix B for details). In other words. Fig. 3 
just shows the bounds of the Markovian region and the non- 
Markovian region. However, the two layers in Fig. 4 corre¬ 
sponding to the different initial states imply the opposite prop¬ 
erties of the acceleration especially for 2 < s < 2.5. Based 
on Fig. 3, the non-Markovianity could effectively signal the 
change of quantum evolution velocity for /3 = ■^, but it fails 



s 


FIG. 3: (Color online). The illustration of the sign of the dephasing 
decay rate 7 ^(t). The green region stands for the positive part of the 
7 '(f) which corresponds to the Markovian region, i.e., the inequality 
Eq. (21) kn < sarctan(f) < {2k -F Ijrr satisfied, and the blank 
region corresponds to the negative part of 'y'{t) which corresponds 
to the non-Markovian region. If f —^ oo, the left critical edge is 
s = 2. 



FIG. 4: (Color online). (including FqI^) versus with the Ohmic 
parameter s and coupling constant rj. The contour corresponds to the 
upper layer. At s ~ 2.5, there is a straight contour, which signals no 
acceleration of quantum evolution for s < 2.5. The actual evolution 
time is f = 3. 


for the lower layer. This shows the strong dependence of the 
initial state. Similarly, the NI-QSF can be given by 


t 


NI _ 
QSL — 


(1 - 

fo |e“^('^i7'(T) I dr 


(23) 


Comparing Eq. (19) with (23), one can immediately find that 
the NI-QSF does not depend on the population of the 
excited state. It is interesting that in this particular model, the 
NI-QSF does not depend on the initial state either. The MF- 
QSF tQgj^ is connected with the NI-QSF by a factor 

2f3y/l — /3^. Thus the NI-QSF versus s and 7 just consistent 
with the MF-QSF for /3 = i.e., the upper layer in Fig. 4. 

Since the opposite prediction on the evolution acceleration has 
been analyzed by the two layers, it is also implied that the NI- 
QSF and the MF-QSF demonstrate contradictory predictions. 


V. DISCUSSIONS AND CONCLUSION 

One may notice that the ML-QSL is zero for some 
special states, such as the ground state |0) for damped J-C 
model, the excited state | 1 ) or ground state | 0 ) for dephasing 
model. However, the ‘distance’ (Bures angle) between the ini¬ 
tial state and some other target state is also zero. So it is not 
difficult to understand the zero evolution time. In addition, the 
NO-QSL is obtained without the von Neumann trace inequal¬ 
ity and operator norm used, so it is obviously tighter than (at 
least, consistent with ) the previous bounds. 

In summary, we have studied the ML-QSL and the NI- 
QSF of the damped J-C model and the dephasing model, and 
found that the importance of the initial-state dependence for 
the quantum evolution acceleration. To some extent, this is 
consistent with the derivation and conclusion the MT-QSL in 
the open system dynamics, which is evaluated with respect to 
the initial state [23]. We also find that the predictions on the 
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evolution acceleration based on different QSLs could produce 
contradictory conclusions. This implies that the QSL could 
deserve us further consideration. We hope that our work could 
deepen the understanding of the quantum speed limit on the 
dynamics of open systems. 


Appendix B: The non-Markovianity of the dephasing model 
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The memory effects (non-Markovianity) for the dephasing 
model associated with the quantum channel capacity ( 5 ($) 
can be defined by [39, 41] 


Appendix A: The non-Markovianity of damped J-C model 

Following Ref. [34], the non-Markovianity measure of the 
damped J-C model we employed is based on the trace distance 
D{^Pi,^P 2 ) = ^tr |4>pi — 4 >P 2 |- The change rate of the 
trace distance is dtD{^pi,^p 2 ). The positive change rate 
stands for the flow of information from the environment back 
to the system. The non-Markovianity of the quantum process 
can be given by 

iV($) = max j dtdtD{^pi,^p 2 ). (24) 

Pi,2(0) JdtD>0 

It is hard to obtain the optimal initial state pair (pi, p 2 ) for 
a general process. However, for the damped J-C model, it 
proves that the excited state | 1 ) and the ground state | 0 ) are 
the optimal state pair, so dtD{^pi, $^ 2 ) can be given by a 
simple expression dtD{pi{t), p 2 {t)) = 2 |g(f)| • \q{t)\' with 
|g(f)f being the time derivative of |(?(f)|. Thus one can eas- 


N = [ dtdtQ{<P), (25) 

JatQ(<s>)>o 





FIG. 5: The non-Markovianity vs. 70 . The parameter A = 15. The 
non-Markovianity measure is normalized to unity. 

ily obtain the relation between the non-Markvianity and the 
parameter 70 , which is plotted in Fig. 5. Comparing with 
Figs. 1 and 2, one can And that the quantum evolution can be 
accelerated in the non-Markovian regime based on NI-QSL. 
However, based on the ML-QSL, the quantum evolution is 
accelerated not only in the non-Markovian regime but also in 
the Markovian regime. 

(a) 



where quantum channel capacity Q{t) = 1 — 772 
with 772 (•) being the binary Shannon entropy. Due to Eq. 
(25), the non-Markovian regime is determined by c)t(5(<f’) > 
0. It is easy to And that dtQ{^) > 0 means 7 '(f) < 0, 
based on which we plot the Markovian and non-Markovian 
regions in Fig. 3. The non-Markovianity versus the param¬ 
eter s and t is plotted in Fig. 6 (a). Comparing with Fig. 
3, one can And that the blank regime in the Fig. 3 corre¬ 
sponds to the non-Markovian regime. In Fig. 6 (b), we plot 
the non-Markovianity versus the Ohmic parameter s and cou¬ 
pling constant p when the evolution time 7 = 3. Comparing 
with the Fig. 4, one can easily And that the quantum evolution 
is accelerated within the non-Markovian regime based on the 
ML-QSL for /? = l/v/2 (the upper layer in Fig. 4). How¬ 
ever, for the lower layer in Fig. 4 = 0.5), the 

quantum evolution can also be accelerated in the Markovian 
regime. This implies the initial-state dependence. 
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